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ABSTRACT 

This paper performs a semi-analytic study of relativistic blast waves in the 
context of gamma-ray bursts (GRBs). Although commonly used in a wide range 
of analytical and numerical studies, the equation of state (EOS) with a constant 
adiabatic index is a poor approximation for relativistic hydrodynamics. Adopting 
a more realistic EOS with a variable adiabatic index, we present a simple form 
of jump conditions for relativistic hydrodynamical shocks. Then we describe in 
detail our technique of modeling a very general class of GRB blast waves with a 
long-lived reverse shock. Our technique admits an arbitrary radial stratification 
of the ejecta and ambient medium. We use two different methods to find dynamics 
of the blast wave: (1) customary pressure balance across the blast wave and (2) 
the "mechanical model" . Using a simple example model, we demonstrate that the 
two methods yield significantly different dynamical evolutions of the blast wave. 
We show that the pressure balance does not satisfy the energy conservation for 
an adiabatic blast wave while the mechanical model does. We also compare two 
sets of afterglow light curves obtained with the two different methods. 

Subject headings: gamma-ray burst: general — hydrodynamics — shock waves 



1. Introduction 

The afterglow emission of a gamma-ray burst (GRB) is believed to be produced by a 
relativistic blast wave (Meszaros & Rees 1997). The relativistic blast wave is driven by an 
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"ejecta" , which is ejected by the central engine of the GRB explosion. As the ejecta interacts 
with a surrounding ambient medium, two (forward and reverse) shock waves develop (e.g., 
Piran 2004). The forward shock (FS) wave sweeps up the ambient medium, and the reverse 
shock (RS) wave propagates through the ejecta. 

As the blast wave has high Lorentz factors 10 2 — 10 3 (e.g., Meszaros 2006), the FS wave 
is highly relativistic and an equation of state (EOS) with a constant adiabatic index 4/3 
may well describe the gas in the FS-shocked region. However, the strength of the RS wave 
varies as the blast wave propagates. In the case of a constant-density ambient medium, the 
RS wave is initially non-relativistic and then transitions to a mildly relativistic or relativistic 
regime (Kobayashi 2000; Sari & Piran 1995). Thus, an EOS with a constant adiabatic index 
is not adequate for the gas in the RS-shocked region; a variable adiabatic index needs to be 
considered to account for change in the gas temperature. 

Although the EOS with a constant adiabatic index has been widely used in analytical 
and numerical studies of relativistic hydrodynamics, it is valid only for the gas of either 
non-relativistic (with the index 5/3) or ultra-relativistic temperature (with the index 4/3). 
The correct EOS for a relativistic ideal gas is formulated in terms of modified Bessel func- 
tions (e.g., Synge 1957), and its equivalent adiabatic index varies from 5/3 to 4/3 as the 
temperature increases. 

As it is not convenient to deal with modified Bessel functions, there has been effort to 
find simpler EOSs that closely reproduce the correct EOS of a relativistic ideal gas. Taub 
(1948) showed that the choice of EOS is not arbitrary and must satisfy a certain inequality 
(Taub's inequality). By taking the equal sign in Taub's inequality, Mignone et al. (2005) 
derived a simple form of EOS that has correct limiting values 5/3 and 4/3. 

The same EOS as in Mignone et al. (2005) was previously introduced by Mathews 
(1971), considering a relativistic "monoenergetic" gas where all particles have the same 
energy. The validity of this EOS was addressed by Blumenthal & Mathews (1976) for the 
cases of both infinite mean free collision times and very short mean free collision times. 
This EOS was also adopted by Meliani et al. (2004) and Mignone & McKinney (2007). 
In particular, Mignone & McKinney (2007) demonstrated in their relativistic numerical 
simulations that use of an EOS with a constant adiabatic index can significantly endanger 
the solution when transitions from cold to hot gas (or vice versa) are present. 

We use the same above EOS in this paper. Following Mathews (1971), we consider 
a relativistic monoenergetic gas and show that it closely reproduces the correct EOS of a 
relativistic ideal gas. Then we use this EOS to find a simple form of jump conditions for 
relativistic hydrodynamical shocks. This simple set of jump conditions applies to shocks of 
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arbitrary strength. 

A short-lived RS was proposed to explain a brief optical flash (Meszaros & Rees 1999; 
Sari & Piran 1999a, 1999b). A dynamical evolution of such a short-lived RS was studied 
analytically by assuming an equality of pressure across the blast wave (Kobayashi 2000; 
Sari & Piran 1995). The RS wave here is short-lived since the ejecta is assumed to have a 
constant Lorentz factor. However, in general, the ejecta is expected to emerge with a range 
of the Lorentz factors. The shells with lower Lorentz factors will gradually "catch up" with 
the blast wave as it decelerates. Thus, the RS wave is long-lived. Such a long-lived RS was 
studied for a power-law ejecta by assuming a constant ratio of the two pressures at the FS 
and RS (Rees & Meszaros 1998). 

In this paper, we present a detailed description of our blast-wave modeling technique 
for even more general class of explosions where the ejecta and the ambient medium have 
an arbitrary radial structure or stratification. More specifically, we study analytically the 
spherical expansion of such a stratified ejecta and find the trajectory of the RS wave through 
the ejecta self-consistently. In order to find a dynamical evolution of the blast wave, we use 
two different methods: (1) customary pressure balance and (2) the "mechanical model" 
(Beloborodov & Uhm 2006). 

Using a simple example model, we demonstrate that, although the customary assump- 
tion of pressure balance for the blast wave yields an estimated evolution, it is not rigorously 
accurate. In particular, the energy conservation is not satisfied for an adiabatic blast wave; 
the total energy is decreased by a factor of 5 in the case of the example model. 

The mechanical model was developed for relativistic blast waves, by relaxing the pressure 
balance (or proportionality) and applying the conservation laws of energy-momentum tensor 
and mass flux on the blast between the FS and RS. Using the same example model, we 
show that the energy conservation is satisfied for the mechanical model. We also show that 
dynamical evolutions found by the two methods differ significantly. Finally, we present the 
afterglow light curves in X-ray and optical bands. We compare two sets of light curves 
corresponding to the two different dynamical evolutions mentioned above. 

In Section |2j we derive a simple set of jump conditions for relativistic hydrodynamical 
shocks. In Section [3j we describe in detail our blast-wave modeling technique. We also 
provide a simple method of evaluating the blast energy, employing a Lagrangian description 
for the blast wave. In Section HJ we review the mechanical model including more detailed 
equations. 
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2. Relativistic shocks 

We consider a shock wave of an arbitrary strength. The preshock medium (cold) is 
denoted by region 1, and the postshock medium (hot) by region 2. The gas moves at right 
angles to the surface of discontinuity (shock front). The rest-mass density p, the energy 
density e (including rest energy), and the pressure p of the gas are defined in the rest frame 
of each region. The energy-momentum tensor T Q/3 for a perfect fluid and the mass flux j a 
are given as 

T afi = (e + p)u a u p + g afi p and j a = pu a , (1) 
where g al3 is the Minkowski metric, and u a is the 4-velocity of the gas. 



2.1. Jump conditions 

A shock is described by three jump conditions that express the continuity of mass, 
energy, and momentum flux densities, respectively, in the shock frame (Landau & Lifshitz 
1959), 

I2P2P2 = IxPiPx, (2) 

TaMea+pa) = 7i7W 2 , (3) 
llPl(e 2 +p 2 )+p 2 = ilftpxc 2 . (4) 

Here subscripts 1 and 2 refer to regions 1 and 2, and c is the speed of light. The Lorentz 
factor 7 of the gas for each region is measured in the rest frame of shock front, and thus 
(3 = (1 — I/7 2 ) 1 / 2 is the gas velocity relative to the shock front. Note that we have set p\ — 
and e\ = pic 2 here since we assume that region 1 is cold. 

We introduce a quantity k in writing an EOS for the relativistic gas in region 2, 

P2 — K 2 ( e 2 ~ P2C 2 ). (5) 

We use the quantity k instead of an adiabatic index, in order to avoid any confusion with 
the Lorentz factors. Note that k 2 is 2/3 for a non-relativistic shocked gas and 1/3 for an 
ultra-relativistic gas. As mentioned in Section [1] we do not use a constant value for k 2 . The 
quantity k 2 here varies between 1/3 and 2/3 as the strength of the shock varies. 

We solve these four Equations fl2J)-([5]) algebraically. First, we define a compression 
ratio a = p 2 j p\ and get the relation 7 2 = a 2 (7| — 1) + 1 by squaring Equation (J2|). Then, 
substituting 7 2 and e 2 (from Equation (jSJ)) into Equations (j3J) and @, we find the expressions 
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for 71,72, and p 2 in terms of k 2 and the compression ratio a, 



2 _ (a + 1) 2 _ (a + l)[a« 2 - (1 + /t 2 )] s 

^2 - G (l-«:|) + (l + / t2 ) 2 ' 71 ~ a(l - K |) + (1 + K2 )2 



2 _ (aK 2 ) 2 - (aK 2 )(2 + K 2 ) 

V2 \P\C ) — 7— r , [I 

(1 + «2j 

where we have used the fact that a(l — k 2 ) + (1 + k 2 ) cannot be zero for a > and 1/3 < 
k 2 < 2/3. Equations © and (J7J) are exact. 

The shock strength may be described by the relative velocity (3\ 2 = — /3 2 )/(l — /3i/3 2 ), 
or by the relative Lorentz factor 7 12 = (1 — /3i 2 ) -1 ^ 2 , which is given by 

712 = (1 - /3i/3 2 )7i72 = 7i72 - [(7? - 1)(7 2 2 - 1)] 1/2 - (8) 

The shock strength 7 12 is then directly calculated using Equation flSJ), 

(^2 - 1) , Q . 
712 = 77— — 7- 9 
(1 + K2) 

Thus, we find a very simple form for the compression ratio a in terms of k 2 and 7i 2 , 

q=^= (1 + " 2)7l2 + 1 . (10) 

Pi «2 

Substituting the compression ratio ( fTUl) into Equations ([6]) and (I?]), we now express 71,72, 
and p 2 in terms of k 2 and 712, 

2 _ (7i2 + !)[(! + ^2)712 -M 2 ri1 x 
71 " (1 - «I)7i2 + (1 + «1) ' 1 J 



2 v ri2 -r ^ / 19 x 

72 (1 - «I)7i2 + (1 + 4) ' 1 J 

P2 = (712-1)[(1 + K2)712 + 1]P1C 2 . (13) 

As we shall see below, n 2 here is in fact a function of 712, and therefore the shock strength 
712 becomes the only free parameter of the shock. Combining Equations (jSJ), (TTOl . and ( ITBl , 
we verify an expected relation, 

e 2 = — p 2 + p 2 c 2 = 712 p 2 c 2 , (14) 

which means that the shock strength 712 is equal to the mean random Lorentz factor of 
particles in the postshock medium (measured in its fluid frame). 

It may be noticed that there are only three independent equations among Equations 
( I10p -(ll4 p . since we had three jump conditions in the beginning. In fact, three of them are 
equivalent to those equations that appear on Blandford & McKee (1976). 
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2.2. Relation between k and mean Lorentz factor 7 

2.2.1. Relativistic ideal gas 



We briefly review a relativistic ideal gas where a particle of mass m and momentum p has 
the energy e(p) = mc 2 [l + (p/mc) 2 ] 1 '' 2 . The Maxwellian momentum distribution function 
f(p) is given as f(jp) = Ap 2 exp [— e(p)/(fesT)], where A is a proportionality constant, fee 
is the Boltzmann constant, and T is the temperature of the gas. Using a normalization 
condition on the number density, n = J °° f(p)dp, the proportionality constant is found to be 
A 



,\ 3 K % u ) , with a modified Bessel function i^" 2 - Here we have defined u = (mc ) / (k B T), 
which basically measures how relativistic the gas is; u ^> 1 corresponds to a non-relativistic 
limit and u < 1 to an ultra-relativistic limit. We evaluate the integrals of pressure and 
energy density of the gas: 



P 



pv(p)f(p)dp = nk B T, 





00 



e(p)f(p)dp 



[rime 



KM | 3 
K 2 (u) u 



(15) 
(16) 



where v(p) is the velocity of a particle of momentum p, and Ki(u) is also a modified Bessel 
function. See Greiner et al. (1995) for an alternative derivation. Equation ( 1T61) gives the 
mean Lorentz factor 7 of particles as 



7 



KM 
KM) 



II, 



The quantity k (defined in Equation (jSJ)) is given as 



(17) 



V 



e — pc A 



u 



KM 

KM 



+ 3 --l 

u 



-1 



18) 



Here the subscript i refers to a relativistic ideal Maxwellian gas. Since both 7 and K{ are 
given in terms of u only, one should be in principle able to express K{ as a function of 7. 
However, it is not easy to do so analytically, as it involves two modified Bessel functions K\ 
and K 2 . This becomes a good motivation of considering the following monoenergetic gas. 



2.2.2. Monoenergetic gas 



Here we consider a monoenergetic gas, where all particles in the gas have the same 
momentum p or the same Lorentz factor 7 (e.g., Mathews 1971). We show below that 
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the behavior of this gas is very close to that of a relativistic ideal gas. The momentum 
distribution function is simply given by a Dirac-delta function, f(p) = n5(p — p), which 
satisfies the normalization condition, n = f(p)dp. We evaluate the integrals of pressure 
and energy density for this gas: 

p = n{mc 2 ) = — , e = n (7772c 2 ). (19) 

37 

Therefore, the corresponding k is found to be 

K m = —^~ 2 =Ul + l), (20) 

e — pc 2 3 \ 7 / 

where the subscript m refers to a monoenergetic gas. Note that K m has the correct limiting 
value 2/3 for a non- relativistic gas and 1/3 for an ultra-relativistic gas. We compare this 
monoenergetic gas with a relativistic ideal gas, by computing n m / Ki numerically as a function 
of 7- The result is shown in Figure [TJ Note that there is only 4.8 % of maximal difference 
at about 7 = 1-6, and two gases are practically identical especially for high Lorentz factors 
above 10. 

Relation (|20p is simple, and thus very useful in dealing with a relativistic gas. In the 
following section, we show that the jump conditions ( |T0|) -(|T4l simplify significantly when the 
gas is treated as monoenergetic. 



2.3. Jump conditions of a monoenergetic gas 

We continue on the problem of a relativistic shock wave for the case of a monoenergetic 
gas. When the gas in the postshock medium (region 2) is treated as monoenergetic, the 
quantity k 2 satisfies 

K2 = \ ( 1 + —), (21) 



3 V 7i2, 

since the mean Lorentz factor in Equation (120]) is equal to the shock strength 7 12 as shown 
in Equation ( Tl4l) . Using relation ( l2~T1) . we rewrite the jump conditions (fT0l -( Tl4"l) in terms of 
712 only: 

ll - (4 1' 1)2 or ft = #2-, (22) 
h 8 7l 2 2 + l W ^2 + 3' 1 J 



2 ^'12 > JX2 mo\ 

72 = 87ITT ° r /?2 = ^' (23) 

P2 = ^(7i 2 -l)PiC 2 , (24) 

a = P2/Pi=47i2, e 2 = 47 2 2 pic 2 . (25) 
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It has become clear that the shock strength 712 is the only free parameter of the shock. 
We emphasize that these simple equations fl22]) - ff25]) are exact for a monoenergetic gas and 
apply to shocks of arbitrary strength (relativistic, mildly relativistic, or non-relativistic). 
This result is also briefly described in Beloborodov & Uhm (2006). 

3. Blast waves 

A central explosion of a GRB ejects a large amount of material with high Lorentz factors 
T e j ~ 10 2 — 10 3 . This ejected flow is called the "ejecta". The ejecta expands and drives a 
forward shock (FS) wave into the external ambient medium. When the ejecta interacts with 
the ambient medium, another shock wave - a reverse shock (RS) - develops and propagates 
through the ejecta. Thus, this standard picture has four regions: (1) external ambient 
medium, (2) shocked external medium, (3) shocked ejecta, and (4) unshocked ejecta (e.g., 
Piran 2004). In Figure |2J we show schematically these four regions. The shocked external 
medium is separated from the shocked ejecta by a contact discontinuity (CD). Two shocked 
regions 2 and 3 between the FS and RS are hot and called the "blast" . 

We assume that the whole blast moves with a common Lorentz factor T. This is rea- 
sonable since internal motions in the blast are subsonic, and hydro dynamical simulations 
confirm that T « const between the FS and RS (Kobayashi & Sari 2000). The Lorentz 
factors Tf, T r , and r e j, denoting for the FS, RS, and ejecta, respectively, are measured in 
the lab. frame. The ambient medium is at rest in the lab. frame. The rest-mass density p, 
energy density e (including rest energy), and pressure p in each region are measured in its 
own rest frame. It is assumed that the ambient medium and ejecta are "cold," having no 
pressure. 

3.1. Radially stratified ejecta 

The explosion ejecta is viewed as a sequence of shells that coast with Lorentz factors r e j. 
Each shell is prescribed an ejection time r. The ejection time plays the role of Lagrangian 
coordinate that labels the shells in the ejecta. Theoretically, the ejecta is expected to emerge 
with a monotonic velocity profile as a result of internal shocks that take place at small radii 
r < 10 16 cm (e.g., Piran 2004). At the end of the internal-shock stage, any two adjacent shells 
no longer collide with each other. Thus, we consider here only a non-increasing function 
of r ej -(r); Pgj(r) = dY^/dr < 0. Note that r e j should remain independent of the lab. 
time t as each shell coasts without colliding with another shell. The corresponding velocity 
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v e j(r) = c (1 — l/T^) 1 / 2 is also non- increasing. The initial density of the ejecta may, however, 
have an arbitrary radial profile. The evolution of such "stratified" ejecta is analytically 
studied here, restricted to the unshocked ejecta. 

We assume spherical symmetry. However, the calculation remains valid even if the 
explosion is driven by a jet with a small opening angle 0j et as long as T e j 3> 0r t ; the jet 
behaves like a portion of spherical ejecta since the edge of the jet is causally disconnected 
from its axis. 



3.1.1. Continuity equation of stratified ejecta 

The 4-velocity u a for a spherically symmetric ejecta is written in spherical polar co- 
ordinates (ct, r, 9, 4>) as u a = T e j (c, v e j, 0, 0), where t indicates the lab. time, and r is 
the radius measured from the center of the explosion. The continuity equation for ejecta is 
simply V a {PejU a ) = 0. Here the ejecta density p e j is measured in the rest frame of ejecta. 
Then the continuity equation reads 



1 d_ 

c dt 



(pejr ej c) 



1 A 



0. 



which becomes 



d_ 

dt 



(p ej r ej ) + ^ 



0. 



(26) 



(27) 



We describe below a simple way of solving Equation ( 1271) that makes use of the Lagrangian 
coordinate r. 

Consider a shell in ejecta that was ejected at time r with velocity v C j(r). The radius r 
of this r-shell at time t is given by 



r{r,t) 



Vej(r) dt' 



,(r) (t-r). 



(2f 



Equation ( 128]) is viewed as a relationship among three coordinates r, t, and r. Any two of 
those can be regarded as two independent variables. In Equation ( 12~T|) . the coordinates r and 
t are two independent variables. Since r e j(r) and v e j(r) are functions of r only, it is useful 
to have r as one of the variables instead of, e.g., time t. Using relation ( )28|) . we eliminate 
time t from Equation ( |27|) . and adopt the coordinates r and r as two independent variables. 



First, we partially differentiate the relation r 

, , , dr 



dr 
dt 







(r) (t — r) with respect to t at fixed r 



+ ^ej(r) 1- 



(29) 
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which yields 



dr 
dt 



(30) 



where vL(t) = dv^/dr. Similarly, we differentiate Equation f l28|) with respect to r at fixed t, 



dr 
Or 



which yields 



Then we find 



dr 
dr 



dr 

, x +^ ej (r) I - — 
Vej{r) V or 



d_ 

dt 

d_ 

dr 



dr 
dt 
dr 
dr 



d_ 

lh 
d_ 
dr~ 



+ 



+ 



dr 
dt 
dr 
dr 



d_ 

dr 
d_ 

dr 



1 



d 



l-r{v' c] /v 2 c] ) dr 
-1 d_ 
Ve] - r{v'Jv e ^ dr 



+ 



d_ 

dr 



(31) 
(32) 

(33) 
(34) 



We substitute Equations ( 1331) and ( J34l) into Equation (1271) . and divide the whole equation 
by Pe^ejVey Then we get 



d 



v ei-r(v' e] /v C j) dr 



in J- 



d_ 

dr 



[lnp e j + ln(r ej w ej )] + - = 0, 



(35) 



where the term ln(r e jt> e j) vanishes since it is a function of r only. Equation f l33|) becomes 



-v' e] /v ei 



+ 



d 



v ej-r(v'Jv ej ) dr 



[lnp ej ]+- = 0. 



Finally, we integrate Equation (|36|) over r at fixed r, 



In 



+ lnp ej +21nr = ln [f(r)} 



(36) 



(37) 



where /(r) is an arbitrary positive function of r. Thus, we find an analytical solution of the 
continuity equation for stratified ejecta, 

n -1 



Pej(r,r) 



f(r) 



r v'~\ 

1 — r-£ 


1 _ f(r) 


\- r - 


v ■ 




c 



C (I*. - 1)3/2 



(38) 



Here we have used the relation v' e Jr) = (c 2 /v c ^)(T' c -/Tl- ). The function /(r) is determined 
below by the initial profile of ejecta near the center of the burst. The factor 1/r 2 represents 
an overall side expansion in 3 dimensional space. The remaining factor inside the brackets 
is responsible for a local spread-out of the ejecta due to its stratification. For 1 dimensional 
plane-parallel symmetry, the solution p e j remains the same as above, except that there is no 
1/r 2 factor. 
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3.1.2. Initial profile of the ejecta 

The r-shell of energy 5E e j(r), ejected at an ejection time r, coasts with its Lorentz factor 
r e j(r). The initial profile of the ejected flow is set by the central engine of the explosion. 
Note that it is completely described by two functions; r ej -(r) and L e j(r) = dE^/dr. Here 
the luminosity L e j(r) is related to the mass flow rate M(r) as L e j = T^Mc 2 . For small radii 
r, the mass flow rate is given by M = (47rr 2 t> e j) (p e j T e j ) . Then the initial profile of ejecta 
density p C j at the burst place is expressed as 

Pej(r,r)= 2 Lej (39) 

The solution (138]) has a limiting form p C j = f(r)/r 2 for small radii r. Thus we determine the 
function /(r), 

/(r) = i^rp = (40) 



Hence the ejecta density is derived as 

LAt) 



47rr 2 v e j r 2 jC 2 



r 



1 - - 



c (r 2 - 1)3/2 



-i 



(41) 



The solution ()4ip is exact; we remark, however, that r^fr) < is assumed in the derivation. 
For given r e j(r) and L e ^(r), the solution fj4T|) allows us to fully understand the subsequent 
evolution of the ejected flow. The solution (jUJ) is presented in Uhm & Beloborodov (2007) 
with no derivation. 



3.2. Jump conditions of the FS and RS 

In Section^ we derived a simple form of jump conditions for shocks of arbitrary strength. 
We apply those jump conditions to the FS and RS of the blast wave, treating the gas in the 
blast as monoenergetic. 

The Lorentz factors 71 and 72 are measured in the rest frame of the FS. As the ambient 
medium is at rest in the lab. frame, we note that 

7i = I7, 72 = (1 - P/3 f )rr f , and 712 = T. (42) 

The relative Lorentz factor 7 12 = T describes the shock strength of the FS. Then the jump 
conditions of the FS read 

2 (4r 2 -l) 2 .4/3 
7? - or A = ^, (43) 
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P2 = ^(r 2 -l) Pl c 2 , /C2 = I^l + Ij > (45) 

p 2 = 4r Pl , e 2 = 4r 2 Pl c 2 . (46) 

Thus, the Lorentz factor Tf and the thermodynamic quantities p 2 , e 2 , P 2, and k 2 immediately 
behind the FS are found in terms of T and p\. Here pi should be evaluated for the ambient 
medium immediately ahead the FS, which we denote by pi(FS). 

The RS is described by the same set of jump conditions when index 1 is replaced by 4 
and index 2 by 3. The Lorentz factors 73 and 74 are measured in the rest frame of the RS, 

73 = (1 - pp r ) rr„ and 7 4 = (1 - PM r ej r r . (47) 

The shock strength of the RS is described by the relative Lorentz factor 743, 

743 = (1 - PM 7473 = (1 - PPei) rr ej . (48) 

Then the jump conditions of the RS read 



2 97 2 3 /3 43 

73 = 87I3TT ° r /?3 = ^' (50) 
4,o^o 1 ( . 1 



P3 = ^(74 2 3 -l)P4C 2 , K3 = ^(l + — ), (51) 

6 6 \ 743/ 

P3 = 4743 p 4 , e 3 = 4743 p 4 c 2 . ( 52 ) 

Equation (147|) yields the Lorentz factor r r , 

r, = (1 - pp 3 ) r 73 = (1 - /3 ej /3 4 ) r ej74 , (53) 

where 73 and 74 are now given in terms of 743 = (1 — Pp e j) rr c j above. Thus, the Lorentz 
factor T r and the thermodynamic quantities ps, e 3 , p 3 , and k 3 immediately behind the RS 
are found in terms of T, r e j, and p 4 . Here T e j and p 4 should be evaluated for the shell 
immediately ahead the RS, which we denote by r e j(RS) and p 4 (RS) = p e j(RS), respectively. 

Let the RS be located at radius r r and sweep up the r r -shell in the ejecta when the 
FS is located at radius 77; the subscripts r and / refer to the RS and FS, respectively. 
When three functions pi(r), r e j(r), and L e j(r) are known, we find then pi(FS) = P i(r/), 
r e j(RS) = r ej -(r r ), and p e j(RS) = p e j(r r , 77); the ejecta density p e j(RS) is given by the solution 
EH) 
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Thus the Lorentz factor F becomes the only free parameter for the blast wave with 
known input functions pi, r e j, and L e j. This is justified since we are given 6 independent 
equations (3 plus 3 jump conditions) for 7 unknowns, which are T, Tf, T r , two independent 
thermodynamic quantities describing the gas behind the FS, and another two for the gas 
behind the RS. 



3.3. Trajectory of the RS through ejecta 



The path of the RS needs to be consistently tracked, as it propagates through the ejecta. 
Consider the RS located at radius r r {t) at time t, sweeping up the r r (t)-shell. Equation 
gives the radius of the r r -shell at time t, which equals r r (t); r r (t) = v C j(r r [ 



we find the velocity v r of the RS, 

dr r 



Using the relation 



dr r 
~dt' 



v e j(r r ) - v ei (r r ) 



1 — 7V 



dr r 



dr r dr r 
dt dt dr r 
we find a differential equation for dr r /dr r , 



dr r 
dr r 



dr r 
dr r 



1 



1 



r? " 2 I 



(t-r r [t\). Then 

(54) 
(55) 

(56) 
(57) 



Equation (|53|) gives the velocity v r in terms of T and r e j(r r ). Thus, Equation ( 1571) allows 
us to follow the trajectory of the RS through ejecta when T is known; for an infinitesimal 
displacement Sr r of the RS, we numerically solve Equation (1571) to find 5r r . 

The Lorentz factor T is determined below by two different methods: (1) customary pres- 
sure balance (see Section [375]) and (2) the mechanical model (see Section H]). We demonstrate 
that the method (1) does not satisfy the energy-conservation law for adiabatic blast waves. 



3.4. Adiabatic blast 

As the blast wave propagates through the ambient medium, the blast grows and the gas 
in the blast evolves hydrodynamically. A simple way of dealing with an adiabatic evolution 
of the blast is described here. Since we treat the gas as monoenergetic, we first study the 
adiabatic process of a monoenergetic gas. 
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3.4-1- Adiabatic process of a monoenergetic gas 

Consider a relativistic monoenergetic gas, which has the pressure p, energy density e, 
mean Lorentz factor 7, volume V, number density n, and particle number N. An adiabatic 
process of the gas is defined by d(eV) = —pdV. When the particle number is conserved, 
we have N = nV = const, or dV = —Ndn/n 2 . Recalling Equation f[T9"j) for p and e of the 
monoenergetic gas, we find 



7 2 - 1 

d(nV •ymc 2 ) = nmc 2 — — — 

which yields 



^ dn 



(5* 



3^ 1 

dj=- dn. (59) 



(7 2 — 1) n 

Here m denotes the particle mass in the gas. We integrate Equation (|59|) to find 

n oc (7 2 - 1) 3/2 and p oc -{f - 1) 5/2 . (60) 

7 

Note that p oc 7 4 is verified for the adiabatic process of an ultra- relativistic gas 7> 1. Using 
Equation ([20]) . i.e., k = ~(1 + I/7), or 7 = 1/(3ac — 1), we re-write relation (1601 as 

/ 1U =Pm{K-)- (61) 

(K- g) 4 

Here we have defined the function p m (K) for the right-hand side. The function p m (K) is 
monotonically decreasing in its valid range, | < k < §. The proportionality constant in 
Equation flBTj) is related to the entropy of the gas, which is conserved for the adiabatic 
process. 



3.4-2. Evolution of adiabatic blast 

We discretize the external ambient medium and ejecta into spherical mass shells 5m, 
and use a Lagrangian description for the blast wave. Each 8m is impulsively heated at some 
point by a shock front (FS or RS), acquiring its initial pressure p and quantity k, which 
are given by the jump conditions (see Section |3~2|) . Using Equation (|6T|) . we can track the 
subsequent adiabatic evolution of each 8m if we know the evolution of its pressure; the initial 
p and k determine the proportionality constant, and we find numerically the new quantity k 
when 5m is at new pressure p. All other thermodynamic quantities of 5m can then be found 
accordingly. For instance, the volume 5V of the mass shell is obtained as 

5V = ^ ^1 5N = ^ - 5iV, (62) 
P 37 p (k-±) 
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using Equation f lT9|) and 7 = 1/(3k — 1). This volume 5V is denned in the rest frame of the 
mass shell. The particle number 5N of the mass shell is calculated when Sm is shocked by 
a shock front, by making use of the particle number flux in Equation ([2]). 



3.4-3. Energy of the blast 



We may then calculate the energy of the blast. When the blast has an instantaneous 
Lorentz factor T, the total energy of the entire blast is evaluated in the lab. frame by 
integrating the 00-component of energy-momentum tensor over the volume of each Sm: 



£ b i ast = ;c[r 2 (e+ P )-p] (^) = E 

{Sm} ^ ' {Sm} 



\ < <)\'+ { r - ^ ) H T 



(63) 



Due to the Lorentz contraction, SV/T is the volume of Sm in the lab. frame. Replacing the 
energy e5V by (7mc 2 ) SN and using Equation (162 p for pSV, we find 



-^biast — 

{Sm} 



ry + Ur- 1 - 
' 3 V r 



7 



7 



(mc 2 <5iV) 



(64) 



We emphasize that the second term here needs to be included in order to correctly express 
the energy of blast. For relativistic blast waves, T 2 3> 1, Equation (164"]) becomes 



E 



Mast — ^ r 
{Sm} 



3 V 7 



(mc 2 5A^) = r 

{<5m} 



^l-/c)(3« + l) 
(3/t- 1) 



(mc 2 5iV) . (65) 



The energy £4 of unshocked ejecta (region 4) is easily found as E4 = L C j(r) dr, 
where r r indicates the location of the RS in the ejecta. The energy of region 1 is negligible 
since the ambient medium is at rest in the lab. frame. Thus, the total energy of the entire 
system is obtained as E tot = E^iast + E4. 



3.5. Customary pressure balance: pf = p r 

A customary approximation assumes a pressure balance across the blast wave; i.e., the 
pressure pf at the FS is equated to the pressure p r at the RS. Equations ( 145]) and ( 151]) give 
the pressures pf = p 2 and p r = p 3 , 

Pf = l - 1) PlC 2 , Pr = ~ ( 74 2 3 - 1) P4C 2 . (66) 



- 16 - 



For relativistic blast waves (r e j ^> 1 and T ^> 1), the pressure balance pf 
the instantaneous T of the blast wave (Beloborodov & Uhm 2006), 

/ \ 1/21 - 1 ' 2 

1 + 2 ^fe) 



Pf = Pr 



p r determines 



(67) 



denoting p 4 = p e j. Recall that p\ 
here; see Section! 



pi(FS), p e j = Pej(RS), and T e j = r e j(RS) should be used 



The solution (1671) indicates that the dynamical evolution of the blast wave is determined 
by purely input parameters pi, p e j, and T e j of regions 1 and 4. The solution T has no 
information on the thermodynamical status of the gas in blast. This observation makes us 
to doubt the validity of the solution f[BT|) . The assumption pf = p r itself is then doubted. 
As we demonstrate below, the solution (1671) in fact does not satisfy the energy conservation 
for adiabatic blast waves. 



3.5.1. Example model 

The initial setup of an explosion is specified by three functions pi, T e j, and L ej , which 
can be arbitrary as long as ^(r) < 0. We consider a simple example model that assumes 

Lej(r) = L = 10 52 erg/s, r cj (r) = 500 - 9r, < r < r b = 50 s. (68) 

The luminosity L e j remains at a constant Lq during the duration t\, of the burst. The Lorentz 
factor r e j decreases linearly from 500 to 50. The total energy Ef, of the burst is simply Ef, = 
LqTi,. The ambient medium is assumed to have a uniform density rii = pi/m p = 1 cm~ 3 . 
Here m p is the proton mass. 

For this example burst, we find the evolution of the blast wave as follows. Suppose 
that, at time t, the RS is located at radius r r (t), the FS is located at radius rj(t), the 
r r (t)-shell passes through the RS, and the blast has the Lorentz factor T(t). We evaluate 
Pi(FS) = pi(r f ), T ej (RS) = r ej (r r ), and p ej (RS) = p cj (r r ,r r ) (Equation flUD), which in 
turn gives 743, 73, T r (Equation (|53l) ). and Tf. For an infinitesimal displacement 5r r of the 
RS, we numerically solve Equation (|57|) to find Sr r . Equation (|28|) . r = f C j(r) {t — r), i.e., 
t = r/f C j(r) + r, gives then the time for the new location of the RS (r r + 5r r in radius and 
r r + 5t t in ejecta). Thus, we find the time interval St for this displacement 5r r . The new 
location of the FS is found by its displacement 5rj during St with the velocity given by its 
Lorentz factor Tf. We evaluate pi(FS), r e j(RS), and p e j(RS) again for the new location. 
Then the solution fl67|) determines the Lorentz factor T of the blast for the new location. 
The result is shown in Figure [31 
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The solution fl67|) also determines the pressure p = pf = p r across the blast. This 
enables us to track the adiabatic evolution of the mass shells in blast (see Section 13. 4. 2p . and 
to find the total energy E tot of the entire system (see Section 13.4.31) . In Figure IU we show 
the resulting total energy E tot . Apparently, the energy conservation is not satisfied; E tot has 
decreased by a factor of 5 by the moment the RS crosses the last shell (r = 50 s) in the 
ejecta. 

3.5.2. What is wrong? 

The spherical expansion of ejecta was completely described by the solution (jUJ. The 
propagation of the RS through ejecta was found self-consistently by Equation ( |5"T|) . The 
conservation laws of energy-momentum tensor and mass flux were explicitly applied across 
both the FS and the RS (the jump conditions). The gas in blast, however, was not required 
to obey those conservation laws; i.e., the pressure balance Pf = p r omits the physics laws 
that should govern the gas in blast. Evidently, this is why the energy conservation for 
the adiabatic blast wave was not satisfied above. Note that one among three independent 
conservation laws was effectively applied to the blast, since we tracked the adiabatic evolution 
of the mass shells in blast in order to find the total energy of the system. 

Adiabatic expansion of the gas in blast implies a pdV work done to the gas itself. This 
work needs to be converted to the kinetic energy of the bulk motion of blast. Clearly, such a 
conversion mechanism is absent from the solution (16?]) as it depends only on input parameters 
of regions 1 and 4. A correct modeling for the blast wave should indicate a mechanism that 
the dynamical variable V is affected by the thermodynamic status of the gas in blast. 

Applying the conservation laws of energy-momentum tensor and mass flux to everywhere 
on the blast, we developed a simple "mechanical model" for the blast wave (Beloborodov 
& Uhm 2006). The mechanical model successfully resolves the energy- violation problem, 
because it replaces the pressure balance pf = p r by the physics laws. We summarize the 
model below including more detailed equations. 

4. Mechanical model for relativistic blast waves 

The gas in the blast wave flows radially with the 4-velocity u a = 7 (1, 0, 0, 0) in spher- 
ical coordinates (ct, r, 9, (ft), where the metric ds 2 = —c 2 dt 2 + dr 2 + r 2 d9 2 + r 2 sin 2 9d(p 2 has 
the determinant g = — r 4 sm 2 9. For any scalar function /, the covariant divergence of fu a 
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is given as (e.g., see Carroll 2004), 

V a {fu a ) = -L^^M^^V) (69) 

where 4 = + is the convective derivative. The rest-mass conservation V a (pu a ) = 
reads then 

1 d f , v. dB 

ssf^^r 1 (71) 

The energy-momentum tensor for a perfect fluid is written as 

T afl = huV + g a y, h = e + p. (72) 

The conservation V^T^ = V ^{hUaU^ + 5£ V) — gives two independent equations (a = 0, 1) 

Id dB 
"9~3I (r 2 h-fu a ) + h'jUa— + d a p = 0, (73) 
r^cdt or 

where Equation (170]) is used. For a = 1, Equation (1731 becomes 

r 2 c dt^ r ^ dr ^ ^ dr ^ ^ 

Instead of V M T M = (a = 0), we use the projection u a V ^TJ 1 = 0. Since u a u a = —1 and 
u a V^u a = 0, the projection becomes 

V,(hu»)=u a V aP , (75) 

which yields 

r 2 cdt c dt dr 

We apply three independent equations (1711 . ( 1741) . and (176]) to the gas between the FS and 
the RS, and make the approximation 

7(t, r) = r(t), <9/3/<9r = 0, r> < r < r f , (77) 

where r r (t) and are the instantaneous radii of the RS and FS, respectively. Then the 
integration of three equations (ITT]) . (174]) . and (1761) over r between r r and rf (at £ = const) 
yields 

(r 2 S T) - T [ Pr (/3 - B r ) + - P)] = 0, (78) 

^ (r 2 HT 2 B) - T 2 /3 [/i r (/3 - B r ) + h f (B f - B)] = p r - p f , (79) 

r[p r (/3-/3 r )+P/(/3/-/3)], ' (80) 
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where £ = pdr, H = J^ f hdr, P = pdr, c(3 r = dr r /dt, and c(3f = drf/dt. Here we 
have used an identity for any function f(t,r), 



rj(t) 



r r {t) 



d 

dt 



f(t,r) 



dr 



d_ 

dt 



r f (t) 



f{t,r)dr 



r r (t) 



\fr(t)W) ~ &(*)} + ff(t){W) - /?(*)}] 



where f r (t) = f(t,r r (t)) and f/(t) = f(t,rf(t)). The relativistic blast is a very thin shell, 
Tf — r r ~ r/r 2 <C r, so we used r / « r r « r when calculating the integrals. 

We simplify Equations (!78 l -( !8~0l) by making use of T 3> 1. The jump conditions at the 
FS give V) = 2T 2 , (3 f - (3 = l/(4r 2 ), e f = 3p f , and h f = e f + p f = Ap f > p f c 2 . The 
convective derivative d/dt = cf3d/dr may be replaced by cd/dr everywhere, and T 2 /3 by T 2 
in Equation (1T9"]) . Then we get 



r 2 dr ^ 

~ r) 

The FS jump conditions give 



Pr ((3-p r )r 2 + - Pf , 

h r {P-P r )T 2 +p r , 

dP 3 

r 2 ^ + (/ ir -p r )(/3-/3 r )r 2 + ^ / . 



(82) 
(83) 
(84) 



Pf = 4r Pi> 



Pf 



r p x c . 



(85) 



Recalling /5 3 = (/3 — (3 r )/(l — (3(3 r ) and the RS jump condition /3 3 = (3^/3, we find 

P43 043 



0-Pr 



r 2 (3-/3/3 43 ) r 2 (3-/343)' 



(86) 



where the second equality is valid for 743 <C T. For the relativistic blast wave, the Lorentz 
factor 743 = (1 — (3(3 e j) TT e j becomes 



743 



1 /r 



+ 



Then we find 



P-Pr 



43 



r 2 - r 2 

ej 



r 2 - r 2 
r 2 . + r 2 ' 



2r 2 (r 2 +2r 2 )' 



The jump conditions at the RS are 



^7) 



Pr = 4743 Pej , e r = 4743 PejC 2 



Pr = g (743 - !) PejC 2 , 



(89) 
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which yield 

ft = 2 {t + Q^ (90) 

4 /r 2 - r 2 \ 

hr = 3 (jf + ^ + lj p CJ c 2 . (91) 

Here Equation (157)1 has been used. 

This leaves four unknowns in Equations (IS2"j) - (IMj) : E, if, P, and T. One more equation 
is required to close the set of coupled differential equations. We propose the following 
approximate relation: 

H -Y,c 2 = 4P. (92) 

As shown in Beloborodov & Uhm (2006), it is accurate in the limits of both an ultra- 
relativistic RS and a non-relativistic RS, and should be a reasonable approximation in an 
intermediate case. 



5. Discussion 

A dynamical evolution of the blast wave is found for the mechanical model as follows. 
For an infinitesimal displacement of the blast, we numerically solve the coupled equations 
(152l)-(l84"l) and ( |92l) of the mechanical model, and find the instantaneous Lorentz factor V and 
integrated quantities H, S, and P. We also solve the differential equation fl5T|) to get the 
RS trajectory through the ejecta. 

For the same example model as used for the customary pressure balance p r = Pf, we 
find the blast-wave evolution. Recall that the example burst is described in Equation (I6"5|) 
and the ambient medium is assumed to have the density rii = pi/m p = 1 cm~ 3 . The result 
found for the mechanical model is shown in Figure |5] (in solid lines); for comparison, the 
solution found for the pressure balance is shown together (in dotted lines). Note that two 
sets of solutions differ significantly; in particular, the blast wave obtained for the mechanical 
model decelerates slower and propagates farther (see Panel c) until the RS arrives at the 
same last shell (r = 50 s) in the ejecta (see Panel a). For the pressure balance, on the other 
hand, the energy loss shown in Figure H] is responsible for the earlier deceleration of the blast 
wave. 

The energy of the blast is easily found for the mechanical model; it is evaluated in the 
lab. frame by integrating the 00-component of energy-momentum tensor over volume, 

Pbiast = I 1 ' [T 2 (e + p)- p] (4irr 2 dr) ~ Aixr 2 (T 2 H - P). (93) 
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We find this blast energy for the same example model above and show the result in Figure [6J 
The total energy E tot of the entire system is precisely conserved; the mechanical model is 
indeed a successful remedy for the energy-violation problem. 

The energy of the blast can alternatively be found by using expression (161)1 . In order to 
track the adiabatic evolution of the mass shells in blast (see Section I3.4.2[) and to find the 
blast energy, we need to know a pressure profile for the mass shells. In case of the mechanical 
model, an instantaneous pressure profile for the blast may be approximated by a quadratic 
function of r (ry < r < rf), which (1) matches two boundary values (p r at r r and pj at rf) 
and (2) satisfies the integrated pressure P. The boundary values are met by a quadratic 
function p(r) = a(r — b) 2 + c with 

b - sC'+'J-s^' < 94 > 

The remaining unknown a is determined such that P = J^ f p(r) dr. An upper bound for a 
is required since p(r) is positive everywhere on the blast; a < {-Jpj + a/p^) 2 /(^/ — r r ) 2 . 

We find this quadratic profile for the dynamical evolution shown in Figure [5] (i.e., the 
solid lines found for the mechanical model) and evaluate expression ( 1611) to find the blast 
energy. The result is shown in Figure [71 The total energy E tot is conserved within about 
5 % for this example. Thus, the quadratic pressure profile should be a reasonably good 
approximation for the mechanical model. 

We also calculate the synchrotron emission from both the FS-shocked and RS-shocked 
regions. We make use of the standard prescription of microphysical parameters (e.g., Piran 
2004): e e (fraction of the shock energy that goes to electron acceleration), e# (magnetic 
parameter), and p (slope of the electron spectrum). We track the synchrotron emissivity 
of all shells on the blast; i.e., a radiative and adiabatic cooling of the electron spectrum 
is tracked for each shell, and an adiabatic evolution of shocked gas is tracked to give an 
evolution of the magnetic field for each shell. The velocity and spherical curvature of the 
shells are also taken into account. A more detailed description will be presented elsewhere. 
The resulting afterglow light curves are shown in Figure El (R band) and Figure [9] (1 keV). 
We find that the two different blast-wave evolutions shown in Figure yield significantly 
different sets of light curves for the same example burst. In particular, the light curves 
obtained with the pressure balance decrease earlier than those with the mechanical model, 
due to an energy loss seen in Figure HI 

A spherically symmetric formulation presented here becomes less accurate when the 
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sideways expansion becomes important at late stages of the blast- wave evolution; i.e., F < 
9^1- In their relativistic hydrodynamical simulations, Meliani et al. (2007) showed that a 2D 
jet-like model decelerates earlier than its ID isotropic counterpart when thermally induced 
expansions lead to significantly high lateral speeds. On the other hand, Zhang & MacFadyen 
(2009) showed that the sideways expansion is a very slow process and previous analytic works 
(Rhoads 1999; Sari et al. 1999) greatly overestimated the rate of the sideways expansion. 

6. Conclusion 

As the blast wave propagates, the strength of the RS wave exhibits a transition from 
non-relativistic to mildly relativistic or relativistic regime (or vice versa). Thus, an EOS 
with a constant adiabatic index is not adequate for the RS-shocked region. We address that 
a more realistic EOS with a variable adiabatic index needs to be used for the gas in the 
RS-shocked region. 

Following Mathews (1971), we consider a relativistic monoenergetic gas and find its 
EOS. We show that there is only 4.8 % of maximal difference in the quantity k (pressure 
divided by internal energy density) when compared to a relativistic ideal gas. Then we 
show that jump conditions of relativistic hydrodynamical shocks simplify significantly for 
the monoenergetic gas (see Section [2]). The simple form of jump conditions presented here is 
exact for a monoenergetic gas and applies to shocks of arbitrary strength (relativistic, mildly 
relativistic, or non-relativistic). We emphasize that its usage is not to be restricted to GRB 
blast waves; it can be applied to other areas of relativistic hydrodynamical shocks. 

Then we present a semi-analytic formulation for relativistic blast waves with a long-lived 
RS. We describe in detail a complete set of tools for finding a dynamical evolution of the 
blast wave for a very general class of explosions. The ambient medium can have an arbitrary 
radial profile, and the explosion ejecta can also be arbitrary as long as r^(r) < 0. We 
provide two different methods of finding dynamics of the blast wave: (1) customary pressure 
balance and (2) the mechanical model (Beloborodov & Uhm 2006). Using a simple example 
model, we show that the pressure balance across the blast wave does not satisfy the energy 
conservation for an adiabatic blast wave; the total energy is decreased by a factor of 5 in the 
case of the example model. 

The mechanical model does not assume a pressure balance or proportionality across the 
blast wave (neither Pf = p r nor Pf/p r = const, is assumed). Instead, it finds the dynamics 
of the blast wave from a set of coupled differential equations that express the conservations 
of energy- momentum tensor and mass flux applied on the blast between the FS and RS. 
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Using the same example model, we show that the energy conservation is satisfied for the 
mechanical model as expected. 

We also show that the two methods yield very different dynamical evolutions of the blast 
wave and, as a result, very different afterglow light curves. We conclude that the customary 
prescription of pressure balance poorly describes the dynamics of the blast wave with a 
long-lived RS and resulting afterglow light curves are inaccurate in a significant manner. 
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Fig. 1. — Ratio K m /Ki as a function of the mean Lorentz factor 7 of gas particles. The 
quantity k is defined in Equation (jSJ), an equation of state (EOS) of a relativistic gas. The 
subscripts i and m refer to a relativistic ideal gas and a monoenergetic gas, respectively. An 
expression for each k (/% or K m ) is derived in Section 12.21 Note that there is only 4.8 % of 
maximal difference between the two at about 7 = 1.6. 
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Fig. 2. — Illustrative diagram of 4 regions in a blast wave. The forward shock (FS) sweeps 
up the external ambient medium (region 1) and the reverse shock (RS) propagates through 
the ejecta (region 4). The shocked ambient medium (region 2) is separated from the shocked 
ejecta (region 3) by a contact discontinuity (CD). Two shocked regions 2 and 3 between the 
FS and RS are "hot" and called the blast. The pre-shock regions 1 and 4 are "cold," having 
no pressure. The entire blast is assumed to have a common Lorentz factor V. 
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10 16 10 17 
r r [cm] 

Fig. 3. — Numerical solution for the blast- wave driven by the example burst specified in 
Equation ([68]); L cj (r) = L = 10 52 erg/s and r ej (r) = 500 - 9r for < r < r 6 = 50 s. The 
ambient medium density is assumed to be n\ = pi/m p = 1 cm -3 . Here m p is the proton 
mass. This solution is found as described in Section 13.5. 1[ using Equation f )BT|) of the pressure 
balance p r = pj. Panel (a) shows the r r -shell passing through the RS at radius r r . Panel 
(b) is the ejecta density n e j(RS) = p C j(RS)/m p of the r r -shell. Panel (c) shows the Lorentz 
factor r e j(RS) of the r r -shell and T of the blast, which yields the relative Lorentz factor 743 
in panel (d). Panel (e) shows that a pressure balance p = p r = pj is assumed across the 
blast. However, this numerical solution does not satisfy the energy-conservation law for the 
adiabatic blast wave; see Figure HI 
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Fig. 4. — Energy £7biast of the blast (dotted line), the energy E 4 of region 4 (dashed line), 
and the total energy E tot of the entire system (solid line; E tot = -EWast + E4) for the blast- 
wave evolution shown in Figure [3j The dynamical evolution in Figure |3] is found by using 
the pressure balance p r = Pf for the example burst specified in Equation (I68I) ; L C j(r) = L = 
10 52 erg/s and r e j(r) = 500 — 9r for < r < T b = 50 s. The ambient medium density is 
assumed to be n\ = pi/m p = 1 cm -3 . The energy E b here is the total energy ejected by the 
burst; E b = Lqu = 5 x 10 53 erg. The energies -Ebiast, E4, and E tot are shown in the units of 
Eb for the r r -shell, i.e., the location of the RS in the ejecta. We precisely track the adiabatic 
evolution of the mass shells in the blast (see Section I3.4.2[) and find the total energy E tot 
of the entire system (see Section r3.4.3l) . However, the resulting total energy is clearly not 
conserved above; it has decreased by a factor of 5 by the moment the RS crosses the last 
shell (r = 50 s) in the ejecta. This demonstrates that the solution (I6T|) derived from the 
pressure balance Pf = p r violates the energy- conservation law significantly for the adiabatic 
blast wave. 
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Fig. 5. — Numerical solutions for the blast-wave driven by the same example burst described 
in Equation (|68|) . The ambient medium density is also the same; ri\ = pi/m p = 1 cm -3 . 
The solid (blue) curves are calculated using the mechanical model (see Section |4j. The 
dotted (red) curves show, for comparison, the solution of Figure [3] (found for the pressure 
balance). Two sets of solutions differ significantly; in particular, the blast wave found for the 
mechanical model decelerates slower and propagates farther (Panel c) until the RS arrives 
at the same last shell (r = 50 s) in the ejecta (Panel a). The solid curves satisfy the 
energy-conservation law for the adiabatic blast wave; see Figure [61 
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Fig. 6. — Energy .Ebiast of the blast (dotted line), the energy E$ of region 4 (dashed line), and 
the total energy E tot of the entire system (solid line; E to t = ^biast + E4) for the blast-wave 
evolution shown in Figure [5] (i.e., the solid blue curves found for the mechanical model). The 
energy E b is the same as in Figure HI The energies -Ebiast, E4, and E tot are shown in the units 
of E b for the r r -shell, i.e., the location of the RS in the ejecta. The blast energy -Ebiast here is 
found as in Equation f )93|) for the mechanical model. The total energy is precisely conserved 
above; thus, the mechanical model successfully resolves the energy-violation problem seen in 
Figure |H 



31 



0.8 



n i i i i r 



J tot 



0) 

C 

CD 



0.6 



0.4 



v£ 4 



"E 



blast 



0.2 







J I I I I I L 



10 20 30 40 50 



r r [s] 



Fig. 7. — Same as in Figure EJ except for an alternative method of finding the blast energy 
-Ebiast- The blast energy here is found by evaluating expression flM|) . while making use of 
an approximate pressure profile of a quadratic function for the blast (see Section 5; 4th 
paragraph). The total energy E tot is conserved within about 5 % above; thus, the quadratic 
pressure profile is a reasonably good approximation for the mechanical model. 
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Fig. 8. — Afterglow light curves in i? band, obtained for the same example burst as in 
Figure |5j FS and RS indicate the emissions from region 2 (FS-shocked) and 3 (RS-shocked), 
respectively. The solid (blue) curves are calculated using the mechanical model, correspond- 
ing to the solid (blue) curves in Figure EJ The dotted (red) curves are obtained using the 
pressure balance, corresponding to the dotted (red) curves in Figure |5j The emission pa- 
rameters are e# = 0.01, e e = 0.1, and p = 2.3. The burst is assumed to be located at a 
cosmological redshift z = 1. The two different blast-wave evolutions shown in Figure |5] yield 
significantly different light curves. The afterglow calculations are terminated when the RS 
arrives at the last shell (r = 50 s) in the ejecta. Thus, from t b s ~ few x 10 3 s, the light 
curves are produced by high latitude emissions. 
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